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1  Software  Description 

The  software  consists  ofDownscaleConcept.exe,  an  executable  file  developed  from  Matlab,  and 
a  GIS  tool  for  running  DownscaleConcept,  DoConGIS,  which  is  used  from  within  ArcGIS. 
DownscaleConcept  was  developed  to  estimate  soil  moisture  at  a  fine  resolution  based  on  coarser, 
spatial  soil  moisture  estimates  by  considering  terrain  characteristics  such  as  slope,  aspect,  and 
contributing  area. 
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2  Software  Installation 


2.1  Software  Requirements 

The  software  was  developed  using  a  Windows-based  computer  and  should  be  run  on  a  computer 
using  either  Windows  XP  or  Windows  7  but  may  be  either  a  32-bit  or  64-bit  architecture. 
Additionally,  the  Matlab  Compiler  Runtime  (MCR)  is  required.  This  software  is  included  in  the 
self-extracting  installation  package.  DownscaleConcept  may  be  run  from  the  Windows  command 
prompt  but  ESRI  ArcGIS  9.3  or  higher  is  required  in  order  to  use  the  DoConGIS  tool.  Two 
versions  ofDownscaleConcept.exe  are  provided:  DownscaleConcept32.exe  is  for  computers 
with  a  32-bit  processor,  and  DownscaleConcept64.exe  is  for  computers  with  a  64-bit  processor. 

2.2  DoConGIS  Installation 

A  self-extracting  *.msi  file  is  provided  for  automatic  installation  of  both  DownscaleConcept.exe 
and  the  DoConGIS  tool.  Double-clicking  on  the  *.msi  fde  will  lead  the  user  through  the 
installation  process.  If  MATLAB  R2010b  or  MCR  (version  7.14)  is  already  installed  on  the  host 
machine,  then  press  Cancel  when  the  following  box  appears.  Otherwise,  this  will  guide  the 
installation  of  MCR  version  7.14. 


When  finished  with  the  installation  process,  the  following  dialog  should  let  you  know  that  the 
installation  was  successful.  Section  3.2  below  describes  how  to  open  the  tool  within  ArcMap. 
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2.3  Optional  Software 

While  not  required  for  use  of  the  downscaling  software,  the  topographic  analysis  freeware 
TauDEM  is  useful  for  calculating  topographic  attributes  and  generating  the  files  needed  as  input 
to  the  downscaling  software.  TauDEM  is  available  at  the  following  URL: 

http :  / /hydrolo  gy.usu.  edu/taudem/taudem5 . 0/ index  .html 

Appendix  B  describes  how  to  acquire  TauDEM  software.  Appendix  C  details  how  to  use 
TauDEM  to  provide  inputs  for  DownscaleConcept.exe. 
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3  Software  Use 


3.1  Inputs 

Use  of  the  DownscaleConcept  software  requires  the  user  to  provide  a  coarse-scale  raster  for  soil 
moisture  as  well  as  fine-scale  rasters  for  elevation,  slope,  specific  contributing  area  (SCA),  and 
aspect  (direction  of  steepest  descent).  Such  terrain  characteristics  can  be  estimated  from  digital 
elevation  models  (DEMs).  Appendix  A  details  the  process  to  obtain  fine-scale  DEMs  from  the 
National  Aeronautics  and  Space  Administration  (NASA),  although  many  other  agencies  provide 
DEM  data.  Additionally,  the  user  must  provide  an  estimate  of  the  annual  average  daily  potential 
evapotranspiration  rate. 

3.1.1  Coarse-scale  moisture  grid 

The  first  input  required  by  DownscaleConcept  is  a  coarse-resolution  soil  moisture  grid.  This 
grid  provides  the  base  moisture  values  that  act  as  spatial  average  moisture  values  for  the 
downscaled  moisture  pattern.  The  moisture  values  in  the  coarse  grid  should  have  units  of  [L3/L3] 
and  the  grid  should  be  based  in  a  projected  coordinate  system  so  that  the  grid  has  horizontal  units 
of  [L],  DownscaleConcept  also  requires  that  the  grid  be  in  a  GeoTIFF  format  (NoData  values 
must  be  a  negative  number). 

3.1.2  Fine-scale  topographic  grids 

DownscaleConcept  requires  fine-scale  grids  for  the  elevation,  slope,  aspect,  and  specific 
contributing  area  (SCA).  The  resolution  of  these  grids  determines  the  resolution  of  the 
downscaled  moisture  estimates  provided  by  DownscaleConcept.  The  extent  of  the  fine-scale 
estimates  is  the  same  as  that  of  the  coarse-scale  moisture  grid.  All  of  these  grids  should  be 
provided  in  a  projected  coordinate  system  so  that  the  grids  have  horizontal  units  of  meters.  The 
DEM  should  have  units  of  [m],  the  slope  grid  should  have  units  of  [m/m],  the  aspect  should  have 
units  of  [degrees  clockwise  from  North],  and  the  SCA  should  have  units  of  [m  / m].  The 
attributes  will  have  the  appropriate  units  when  a  DEM  with  appropriate  projected  coordinates  is 
analyzed  for  slope  and  SCA  by  the  TauDEM  software  and  aspect  is  calculated  by  the  ArcGIS 
Aspect  function  (Spatial  Analyst  toolbox,  Surface  tools).  DownscaleConcept  requires  the  grids  to 
be  in  the  GeoTIFF  format  (NoData  values  must  be  a  negative  number),  which  is  also  the  file 
format  used  as  input  and  provided  as  output  by  the  TauDEM  software.  Instructions  for 
calculation  of  terrain  characteristics  are  found  below  in  Appendix  C. 

In  order  to  minimize  the  impacts  of  cross-boundary  features  (e.g.  watersheds)  on  the  topographic 
attributes  calculated  for  the  area  of  interest,  the  initial  DEM  used  to  calculate  the  attributes 
should  extend  some  distance  beyond  the  area  of  interest  on  all  sides  during  the  DEM  topographic 
analysis.  The  DEM  can  then  be  clipped  to  the  area  of  interest  after  the  analysis.  The  user  may 
need  to  connect  (“mosaic”)  multiple  smaller  DEMs  together  to  obtain  a  single  DEM  extending 
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beyond  the  area  of  interest.  Note,  however,  that  DEM  extents  that  are  too  large  and/or 
resolutions  that  are  too  fine  may  lead  to  computer  memory  problems. 

3.1.3  Potential  evapotranspiration 

The  final  input  required  for  DownscaleConcept  is  an  estimate  of  the  annual  average  daily 
potential  evapotranspiration  (PET)  rate  in  [mm/day]  for  the  estimate  area.  Only  one  PET 
estimate  is  required  for  the  entire  estimate  area. 

3.2  Running  DownscaleConcept 

Instructions  for  using  DownscaleConcept  both  through  the  DoConGIS  tool  and  from  the 
command  prompt  are  provided  below.  The  screen  capture  figures  in  the  instructions  are  the  same 
as  those  the  user  should  see  when  using  the  example  files  provided  in  the  “C:\Program 
Files\CSU\DownscaleConceptGIS\Example”  folder  created  during  software  installation. 

3.2.1  DoConGIS  tool 

1.  Open  ArcMap  after  installing  DownscaleConceptGIS  (refer  to  the  installation 
procedure  asdescribed  in  Section  2.2 above). 

2.  On  the  menu  bar  at  the  top,  click  Tools->Customize...  On  the  Toolbars  tab,  check 
"DownscaleConceptTools" asshown  in  the  figure  below.  ClickClose. 


Toolbars  |  Commands  |  Options  | 
Toolbars: 


□  Data  Frame  Tools 

□  Dimensioning 

G  Distributed  Geodatabase 


□  DownscaleConcept  Tools 


0  Draw 

□  ERA  Aggregator 
0  Editor 

□  Effects 

□  GPS 

□  Geo  Dataset  Exchange  Tools  3 

□  Geocoding 

□  Geodatabase  History 

□  Geometric  Network  Editing 


New1... 


Rename.. 


tl 


Delete 


Reset... 


Keyboard. 


Add  from  file... 


□ose 
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3.  The  following  toolbar  should  now  be  visible. 


4.  Clickon  the  DoCon  icon  and  the  following  dialog  should  appear. 


5.  Specify  the  estimated  annual  average  daily  potential  evapotranspiration 
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6.  Provide  the  77 FF  file  names  for  the  coarse-scale  moisture  grid  and  the  fine-scale 
topographic  attributes  (which  can  be  calculated  from  TauDEM  according  to 
instructions  found  in  Appendix  C)  to  the  dialog  box  either  by  dragging  and  dropping 
raster  layers  into  the  drop-down  boxes,  expanding  the  drop-down  menu  and 

selecting  a  file  path,  pressing  the  J  icon  and  browsing  for  the  files,  or  placing  the 
file  path  directly  into  the  drop-down  box. 
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7.  The  completed  dialog  box  should  look  similar  to  the  one  below  when  using  the  input 
files  from  the  Example  subdirectory  of  the  installation  directory. 


11 


0  I 

8.  Click  the  Save  (ml)  button  on  the  dialog  to  save  your  soil  moisture  model  file  (using 
*.sme  extension).  This  saves  a  text  file  containing  Hie  paths  to  all  input  files  for 
Downscale  Concept aswell asmodel parametersand  running  options 


Save  As 

^l  C:\Program  Files  [x86)\CSU\Do wn seal eConceptGISVExam pie  ▼ 

*t| 

II  Scorch  Example  p  | 

Organize  ■ 

w  Newfolder 

i==  -  # 

I C  Desktop 
l  Downloads 
Recent  Places 

a  Libraries 
-  Documents 
Music 

[^1  Pictures 
0  Videos 

Computer 
i  pi  Local  Disk  (C:) 
rn  nri-rQ  m.’i 


Name 


File  name:  Demo_Dry.sme 


Save  as  type:  File*.sme 


*  Hide  Folders 


Date  modified  Type 

No  items  match  your  search. 


Size 


Save 


Cancel 


9.  Opening  a  SME  file  in  a  text  editor  will  return  something  like  the  following  figure, 
where  the  lines  contain  the  version  number,  file  paths  to  input  files  and  model 
para  meters  or  options  respectively. 


'ffl  I  □  -*>  C 


Demo_Dry.sme  -  WordPad 


0 


Home 


□ 

Paste 


JC  Cut 

4  c°py 


Clipboard 


View 


Courier  New 

B  I  U 


-  11 
K,  X1 


Font, 


\£’\A’\  d \E 


Paragraph 


Picture  Paint  Date  and 
T  drawing  time 

Insert: 


a 

Insert 

object 


© 

£6  Find 
Replace 
[□]  Select  all 
Editing 


[2.0 

C:\Program  Files 
C:\Program  Files 
C:\Program  Files 
C:\Program  Files 
C:\Prograra  Files 
C:\Program  Files 
C:\Program  Files 
C:\Program  Files 
C:\Program  Files 


(xB  6) \C 3U\Downs c a 1 eC one eptGI 3 \ Downs c a 1 eC oncept 6 4 . exe 
(xB  6)  \C 3U\ Downs caleConceptGI3\Pa rams  .  txt 

(xB6>  \CSU\Down3caleConceptGIS\Exaiaple\Coar3eSM_Dry_Prj  .tif 

fxfl  6) \CSU\DownscaleConceptGIS\Example\DEM. tif 

(xB  6) \CSU\Down3caleConcepfcGIS\Example\ Slope .tif 

(xB  6)  \C3U\DownscaleConceptGI3\Example\ContribArea  .tif 

(xB  6)  \C  3U\Downs  c  aleC  onceptGI  S \  Ex  amp  1  e \Aspe  ct  .tif 

(xB  6)  \CSU\DownscaleCoiiceptGIS\Exaniple\PSRI_Dry .  tif 

fxB  6)  \CSU\Down3caleConceptGIS\Exainple\FineSM_Dry.  tif 


1.9128  944.797  5.8419  0.5  0 . 3651  4.134-0.3743 

0.01  0.0135  0.0166666666666667  4.8652  0 


100%  ©  ' - 0 -  © 
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10.  Under  the  Tools  menu  at  the  top  of  the  dialog ,  other  model  parameters  can  be 
viewed  and  adjusted. 

a.  Tools->Units Description...  describesthe  unitsused  foranalysis. 


Units  Description 


0 


“GeoTIFF  Files 
Coarse  Soil  Moisture 
Digital  hlevation  Model  [m] 

Slope  [m/m] 

Aspect  [degrees  clockwise  from  North] 
Contributing  Area  [m*2>!m] 


Other  Parameters 

Potential  Eva  potra  ns  pi  ration  (ET)  [mnVday] 

Vertical,  saturated  hydraulic  conductivity,  ksv  [mnVday] 

Florizontal,  saturated  hydraulic  conductivity,  ksh  [mm/day] 

Exponent  parameter  for  lateral  flow,  gamma  [dimensionless] 

Soil  porosity,  phi  [dimensionless] 

Depth  of  soil  for  location  with  no  curvature,  delta  [mm] 

Vegetation  exponent  parameter  for  ET,  beta  [dimensionless] 
hypothetical  curvature  for  which  location  would  have  no  soil,  kO[mA-1 
Minimum  soil  depth  [m] 

Minimum  PSRI  value  [dimensionless] 

Minimum  s  lope  va  I  ue  [m/m] 

~  OK  | 


b.  Tools->Model  Parameters...  allows  the  user  to  view  or  specify  the  model 
parameters  that  are  conceptually  linked  to  physical  processes  affecting  soil 
moisture  levels. 
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c.  Tools->Atthbute  Bounds...  allows  the  user  to  specify  minimum  values  for  the 
Potential  Solar  Radiation  Index  (PSRI),  which  is  calculated  within 
DownscaleConcept,  and  the  slope.  These  attribute  values  are  found  in  the 
denominators  of  calculations  within  DownscaleConcept  so  their  minimum 
values  must  be  greater  than  zero.  A  nominal  minimum  slope  value  is  0.5  / 
cellsize.  In  this  case  the  cell  size  Is 30  m. 


d.  Tools->Set  File  Paths...  allows  the  user  to  specify  the  paths  to  the 
DownscaleConcept  executable  file  and  the  parameter  file  path. 
DownscaleConceptGIS  writes  the  parameter  file  prior  to  running 
DownscaleConcept.exe. 
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e.  Another  option  within  DownscaleConcept  is  to  truncate  downscaled  soil 
moisture  estimations  at  soil  porosity.  This  can  be  done  by  checking  Tools- 
>Truncate  SM  at  Porosity. 


f.  Hovering  over  any  of  the  model  parameters  in  the  above  dialogs  will  prompt 
a  balloon  to  pop  up,  describing  the  parameter. 
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g.  Clicking  About  brings  up  a  description  of  the  tool. 


11.  Save  the  *.sme  file  again.  Then,  click  "Run  Downsc  aleC  oncept"  at  the  bottom  of  the 
dialog  box  to  calculate  the  downscaled  soil  moisture  estimates 
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12.  A  console  window  will  display  while  Downsc  a  leC  one  ept  is  running: 


13.  When  the  tool  finishes  calculating,  the  output  should  look  like  similar  to  the  following 
figure,  although  the  console  window  will  close  immediately  afterfinishing: 


14.  If  the  output  TIFF  files  are  not  automatically  loaded  into  ArcGIS,  then  either  an  error 
occurred  in  DownscaleConcept  or  in  pulling  output  into  ArcGIS.  One  solution  is  to 
make  sure  all  in  put  TIFF  files  are  in  GeoTIFF  format,  exported  from  TauDEM  orArcMap. 
All  fine-scale  topographic  rasters  must  be  exactly  the  same  size,  containing  the  same 
amount  of  rows  and  columns  with  the  same  cellsize  and  spatial  reference.  Also, 
ensure  that  none  of  the  input  files  found  temporary  directory  are  loaded  within 
ArcGIS  (Temp  files  are  copied  to  the  location  shown  in  Tools->Set  File  Paths...,  and 
then  DownscaleConcept.exe  executes  calculations  on  the  temporary  files,  and  then 
temporary  files  a  re  deleted.  This  process  cannot  be  done  when  these  temporary  files 
are  loaded  within  ArcGIS  because  ArcGIS  places  a  lock  on  any  file  open  within  its 
framework). 
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3.2.2  Command  line  version 


Two  executable  files,  “DownscaleConcept32.exe”  and  “DownscaleConcept64.exe”,  are 
distributed  with  DownscaleConceptGIS  and  can  be  found  in  the  “Program  Files\CSU\ 
DownscaleConceptGIS”  folder  on  the  host  machine.  The  appropriate  executable  file  for  the 
central  processing  unit  of  the  host  machine  may  be  run  from  the  Windows  Command  Line 
Prompt  or  by  double-clicking  on  a  batch  file.  To  run  DownscaleConcept  from  the  command  line, 
the  user  must  provide  seven  file  names,  representing  the  model  parameter  file,  four  input 
GeoTIFF  files  and  two  output  GeoTIFF  files.  The  file  names  must  be  provided  in  the  following 
order:  model  parameters,  coarse  soil  moisture,  elevation,  slope,  specific  contributing  area,  aspect, 
fine-scale  PSRI,  and  fine-scale  moisture  estimate. 

3.2.2. 1  Input  file  format 

When  running  DownscaleConcept  from  the  command  line  the  user  must  first  specify  the  name  of 
the  file  containing  the  model  parameters.  This  file  is  a  standard  text  file  and  the  default 
parameter  values  are  provided  in  the  file  named  “Params.txt”  in  the  installation  directory.  The 
format  of  the  file  is  shown  in  the  below  screen  capture  (A  “%”  sign  at  the  beginning  of  a  line 
indicates  a  comment,  and  values  are  delimited  by  white  space,  either  spaces  or  tabs): 


- 

j|  Params.txt  -  Notepad  [  □  |r4s] 

File  Edit  Format  View  Help 

%Parameter  file 
%  1st  Row: 

%  ksv:  vertical  saturated  conductivity  [mm/day] 

%  ksh:  lateral  saturated  conductivity  [mm/day] 

%  gamma:  exponent  parameter  for  lateral  flow 

%  phi:  soil  porosity 

%  dO:  depth  of  soil  for  location  with  no  curvature,  [tit] 

%  beta:  vegetation  exponent  parameter  for  evapotranspi rati  on 

%  kO:  hypothetical  curvature  for  which  location  would  have 

%  no  soil,  [imA-i] 

% 

%  2nd  row: 

%  mrinDepth:  Minimum  soil  depth 

%  minPSRi:  Minimum  PSRI  value 

%  minslope:  Minimum  slope  value 

%  pet:  Estimated  annual  daily  average  Potential  ET  value 

%  for  entire  region  [mm/day] 

%  truncateSM:  0  if  not  truncating;  1  if  user  desires  to  truncate 

%  soil  moisture  values 

%  xl:  unassigned  value  (must  have  a  number  here  though) 

%  x2:  unassigned  value  (must  have  a  number  here  though) 

% 

1.9128  944.797  5.8419  0.5  0.9651  4.134  -0.3743 
0.01  0.0135  0.0166666666666667  4.8652  1  0  0 

4  ► 


The  remaining  file  names  provide  locations  for  the  coarse-scale  soil  moisture  grid,  the  fine-scale 
topographic  grids  and  the  output  files.  The  requirements  for  the  coarse-scale  moisture  and 
topographic  attribute  input  files  are  discussed  in  Section  3.1  above. 

3.2. 2. 2  Command  line  usage 

DownscaleConcept  can  also  be  run  from  a  batch  file  and  an  example  batch  file 
(“RunDownscaleConcept.bat”)  is  included  in  the  installation  directory.  Batch  files  are  simply 
text  files  with  a  list  of  command  line  actions  to  perform.  Command  line  prompts  must  follow  the 
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pattern  below  (brackets  ‘[’  and  ‘]’  indicate  optional  parameters,  if  running  on  a  32-bit  processor, 
replace  DownscaleConcept64.exe  with  DownscaleConcept32.exe): 

[<path\>]DownscaleConcept64.exe  “[<path\>]<Params  text  filename>” 

“[<path\>]<Elevation  GeoTIFF  filename>” 

“[<path\>]<CoarseSoilMoisture  GeoTIFF  filename>” 
“[<path\>]<Slope  GeoTIFF  filename>” 
“[<path\>]<ContributingArea  GeoTIFF  filename>” 
“[<path\>]<Aspect  GeoTIFF  filename>” 
“[<path\>]<OutputPSRI  GeoTIFF  filename>” 
“[<path\>]<OutputSoilMoisture  GeoTIFF  filename>” 

Different  input  file  paths  must  be  separated  with  a  space  at  the  command  line.  If  a  file  path 
contains  a  space  within  its  name  then  the  file  path  must  be  enclosed  in  double  quotation  marks 
(e.g.,  “File  Path  Containing  Spaces.txt”).  For  an  example  of  a  command  line  prompt  that  runs 
DownscaleConcept,  open  “RunDownscaleConcept.bat”  with  any  text  editor.  The  following 
screen  capture  shows  the  command  prompt  window  when  running  the  command  line  version 
using  the  files  in  the  Example  subdirectory.  The  command  prompt  below  shows  a  general  form 
when  the  current  directory  contains  “DownscaleConcept[64  or  32J.exe.”  Quotations  are  only 
required  around  file  paths  containing  spaces. 
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3.3  Outputs 

Outputs  of  DownscaleConcept,  whether  run  through  DoConGIS  or  from  the  command  line,  are 
GeoTIFF  files  for  the  fine-scale  soil  moisture  and  the  potential  solar  radiation  index  (PSRI) 
pattern.  The  fine-scale  soil  moisture  estimate  provides  soil  moisture  estimates  for  the  specified 
date  at  the  spatial  resolution  of  the  provided  topographic  attributes.  The  estimate  has  the  same 
units  as  that  of  the  coarse  moisture  pattern,  [L  /L  ].  The  extent  of  the  estimate  is  that  of  the 
intersection  of  the  coarse-scale  moisture  and  the  topographic  attributes  provided.  A  second 
output  from  DownscaleConcept  is  the  PSRI  pattern  with  the  same  extent  and  resolution  as  the 
other  topographic  attributes.  PSRI  is  calculated  and  used  by  DownscaleConcept  to  generate  the 
fine-scale  moisture  pattern  and  is  provided  as  output  for  reference.  The  PSRI  represents  the 
potential  insolation  available  on  a  sloped  surface  relative  to  that  of  a  horizontal  surface  at  the 
same  location  on  the  winter  solstice  date. 

3.3.1  Viewing  Outputs  in  ArcGIS 

DoConGIS  automatically  loads  output  files  into  ArcGIS  after  running  DownscaleConcept  when 
Tools->Auto  Load  Output  is  checked.  If  for  some  reason  output  files  do  not  load,  viewing  the 
outputs  will  require  some  post-processing.  Since  GeoTIFF  files  do  not  store  a  tag  for  a  NoData 
value  (a  GeoTIFF  file  may  look  as  if  it  contains  a  NoData  value  within  ArcGIS,  but  NoData 
values  are  actually  stored  within  *.aux  files),  a  filtering  process  is  required  for  ArcGIS  to  read 
the  output  rasters.  Instructions  below  discuss  how  to  view  the  data  in  ArcGIS. 

1.  In  ArcMap,  right-click  on  the  data  frame  in  the  table  of  contents  and  select  Add 
Data... 
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2.  Browse  to  the  folder  containing  your  outputs  and  add  the  outputs  to  the  map.  Select 
Yes  when  ask  to  build  pyramids. 
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3.  ArcGIS  will  automatically  assign  large  lower  and  upper  bounds  to  the  outputs.  To  fix 
this,  right-click  on  the  layer  in  the  table  of  contents  to  open  its  Properties  window  and 
select  the  Symbology  tab.  Select  Classified  in  the  field  on  the  left.  Click  Yes  when 
asked  to  compute  a  histogram. 
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4.  PressOK,  and  the  output  may  look  bad  initially...  like  the  below  figure.  Ihisisbecause 
of  the  unrecognized  NoData  value. 


Untitled  -  ArcMap  -  Ardnfo 
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D  E)| 

QQ  “  I 

g  W  #  ¥  4  CF  B 

x®v  & 

Spatial  Analyst  ^ 

Layer:  |  PSRI_Wet2.tif 
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M  lit 

- '*J 

P  Layers 
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■  -3 .402823 e+ 038 
□  -3.402823e+038  - 1.611108 
3  0  PSRI_Wet2.tif 


gfr  ArcToolbox 
g  ^  3D  Analyst  Tools 
g  ^  Analysis  Tools 
Arc  Hydro  Tools 
p  ^  Cartography  Tools 
g  ^  Conversion  Tools 
p  ^  Coverage  Tools 
g  ^  Data  Interoperability  Tools 
El  ^  Data  Management  Tools 
Geocoding  Tools 
g]  ^  Geostatistical  Analyst  Tools 
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5.  DownscaleConcept  nominally  exports  NoData  in  the  GeoTIFF  file  as  a  negative 
number.  In  order  for  ArcG IS  to  recognize  this  number,  a  new  grid  will  be  calculated 
using  a  condition.  There  are  many  ways  to  do  this  in  ArcG  IS,  but  a  simple  way  is  to 
open  ArcToolbox->Spatial  Analyst  Tools->Conditional->Con.  Place  the  output  fine- 
scale  soil  moisture  grid  into  the  first  and  third  fields  and  place  “ Value  >=  0"  in  the 
second  field  like  so: 
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6.  The  following  raster  will  automatically  be  added  to  your  map  session,  which  looks 
much  better. 
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Appendix  A  Acquiring  ASTER  elevation  data 


A  digital  elevation  model  (DEM)  is  used  to  calculate  the  topographic  attributes  required  by 
DownscaleConcept.  ASTER  global  DEM  data  is  available  for  most  parts  of  the  world  and  can 
be  downloaded  from  the  internet  at  the  URL: 


http://www.echo.nasa.gov/reference/astergdem  tutorial.htm. 

For  convenience,  directions  for  downloading  ASTER  data  have  also  been  transcribed  here: 

1.  Go  to  https//y\d^ecjTojTasa_^py^apjl  and  login  to  your  account  (you  must  be  a 
registered  user). 

2.  Select  the  following  options  in  the  search  Page: 

a.  From  Land  Discipline:  Select  [ASTER]  (This  will  populate  the  Category  box 
definitions) 

b.  From  Category  box:  Select  [ASTER  Global  Digital  Elevation  Model  V001] 


Choose  Keywords  for  One  or  More 
Categories  1 

Pick  a  dtsopitmAoptc  [tor  Arwospherisi.  TRAIMJL  than  tom  ffte  tist  of  data  sals. 

For  rrwltspfe  iiqptos  ctaww  one  (qprc  £  data  s ets.  toeo  fhe  nerf  tepw  £  date  seta 
Td  more  ihan  One  date  gef,.  use  CM-cfrcJr  hr  PCs;  App!e-ct(Ck  for  fite&ntoSh. 


Text  Search: 

GoJ 


Help 
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kSTES.  li a  Ete  cassis:  at  ce-d  "Jap^ectised  rc.i- si-ameci 


View  Data  Set  Definitior-  | 

Atmosphere: 


c  Aqua  AflRS/AMSU-A'HSB 

C 

MISR 
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© 
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© 
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0  CERES/ Aqua 

© 
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0  CERES/Tena 

© 

MOPITT 

C  CERES/TFmi 

o 

TOAA,  AMS  U- A/M  SO 

O  ¥\emn  Situ 

© 

OMl/Aura 

0  PSSP 

© 

SAGE 

C  GLAS/ICES-al 

0 

SSMi'l 

O  GOES 

o 

TES/Aura 

ft  HIRDLS/Aura 

o 

TOMS 

ft  HRJR/MRfiVIWIR 

© 

TOVS 

C  UGHTMING 

© 

TRMM 

ft  LIMS 

o 

UARS 

ft  MERRA 
E*y  Discipline 


Choo  se  Data  Set  Keywords  | 
Cryo  sphere:  Land: 

Oceans: 

Solan'Other: 

C  AMSR/AMSR-E  C  AMSR/AMSR-E 

0  AMSR/AMSR-E 

ft  ACRIM 

C  GLASflCESat 

ft  ASTER 

ft  GLASnCESflt 

C  Field/lln  Silu 

ft  MODIS/Aqua 
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O  Landsat  1-5 

0  MOOliT&fte 
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c  sar 

ft  Land-sat  7 
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ft  MODIS/Aqua 
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"By  Discipline"  not  responding?  Use  [ha  non-savasenm  yers"&rc 

O  By  Citegori^sWttJlbutis  Try  Me’  O  By  Keywords 


3.  Choose  search  area  parameters  location  or  region 

4.  Scroll  to  the  bottom  of  the  page  and  select  “START SEARC  H" 

a.  Note  that  you  can  only  search  for  1000  granules  at  a  time.  If  there  are  more 
granules  that  match  your  query,  you  will  be  first  presented  with  the  dataset 
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results  to  your  query.  From  this  page,  check  the  box  next  to  the  ASTER  GDEM 
dataset  and  press  the  "List  data  granule s"  button  to  view  your  result  granule 
listing. 

b.  To  find  more  than  the  1000  granules  returned,  it  is  recommended  that  you 
refine  your  search  parameters 

5.  Proceed  to  "Step  1"  and  select  your  desired  granules  Select  an  area  greaterthan 
your  area  of  interest  because  edge  effects  can  be  significant  while  calculating 
contributing  area.  The  following  figure  displays  the  area  of  interest  compared  to  the 
downloaded  DEM  coverage.  Notice  that  four  DEM  s  were  downloaded  to  coverthe 
area  of  interest. 


28 


6.  Pressthe  ", Add  selectionsto  cart”  button  to  create  your  order  contents. 

7.  Accept  the  Data  Quality  Summa  ry  Disc  la  ime  r. 

8.  Choose  ordering  option  for  data  retrieval.  WIST  will  report  that  the  granules  are  not 
orderable  if  you  have  not  logged  in  as  a  registered  user. 

9.  Check  appropriate  policy  acknowledgements  (required  fields) 
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10.  If  desired,  select  a  different  option  to  applying  the  order  options  to  all  granules  or 
saving  as  order  preference. 

11.  Proceed  and  select  "Next  Step".  The  data  granules  will  now  be  accessible  in  your 
shopping  cart. 

12.  Select  "Go  to  step  2:  OrderForm" 

13.  You  may  be  asked  to  fill  out  address  location  and  billing  information  for  verification. 
No  charge  will  be  associated  with  GDEM  data  archives 

14.  Select  "Submit  Order  Now!"  to  processyour  order. 

15.  Wait  foremail  confirmation  to  retrieve  data  via  FTP. 

16.  Follow  email  directions  to  retrieve  data. 
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Appendix  B  Acquiring  TauDEM  5  software 

The  TauDEM  software  is  free  software  produced  by  Utah  State  University  that  is  useful  for 
calculating  the  topographic  attributes  required  by  DownscaleConcept.  TauDEM  can  be 
downloaded  from  the  following  URL: 

http://hvdrology.usu.edU/taudem/taudem5.0/index.html. 

That  website  also  provides  instructions  for  installing  and  using  the  TauDEM  software.  Note  that 
TauDEM  also  requires  installation  of  a  message  passing  interface  software,  MPICH2,  and  the 
above-referenced  website  also  provides  instructions  for  obtaining  and  installing  that  software  as 
well.  The  TauDEM  instructions  about  adding  the  MPICH  software  to  the  user’s  search  path 
were  unclear.  Therefore,  additional  instructions  for  this  part  of  the  installation  are  found  below: 

1.  Click  Sta rt->Rig ht-c lie k  Computer->Properties->Advanced  System  settings- 

>Advanced  tab->Environment  Variables...  Navigate  within  the  "System  variables"  to 
find  "Path"  in  the  “Variable"  field  asshown  in  the  screenshot  below. 


User  variables  for  hambre 


Variable 

Value 

TEMP 

c:\temp 

TMP 

c:\temp 

New...  ]  [  Edit...  ]  |  Delete 


System  variables 


Variable 

Value 

- 

□S 

Windows_NT 

i — i 

IPath 

C :  \Windows^ystem32;C :  \Windows;C :  \. . . 

lu 

PATH  EDCT 

.  COM; .  EXE; .  BAT; .  CMD; .  VBS; .  VBEj .  JS; . . . . 

PRDCESSOR A. . . 

AMD64 

New...  |  [  Edit .  ]  |  Delete 


OK  |  |  Cancel 


2.  Click  Edit...  below  the  System  variables.  Add  ";C:\  Program  FileS\  MPICH2\  bin"  to  the 
end  of  the  values  in  the  "Variable  value:"  field. 
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If  ArcGIS  was  running  during  the  installation  of  TauDEM,  close  and  reopen  it  before  using 
TauDEM. 
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Appendix  C  Terrain  Characteristics  -  Calculations 

DownscaleConcept  requires  terrain  information  such  as  contributing  area,  slope  and  aspect  to 
downscale  a  coarse  soil  moisture  estimate.  In  order  to  determine  topographical  characteristics, 
DEM  must  be  downloaded,  manipulated,  and  utilized  in  terrain  characteristic  calculations.  A 
sample  tutorial  for  processing  data  through  TauDEM  is  found  in  the  sections  below. 


C.l  Manipulate  DEM  data  for  placement  into  TauDEM 

Prior  to  calculating  terrain  characteristics,  a  few  pre-processing  steps  should  be  considered.  In 
order  to  minimize  the  impacts  of  underestimating  contributing  area  calculated  for  the  area  of 
interest  (e.g.,  having  a  watershed  much  larger  than  the  downloaded  DEMs),  it  is  suggested  that 
the  initial  DEM  extend  beyond  the  area  of  interest  during  the  DEM  topographic  analysis.  The 
DEM  can  then  be  clipped  to  the  area  of  interest  after  the  analysis.  The  user  may  need  to  connect 
(“mosaic”)  multiple  smaller  DEMs  together  to  obtain  a  single  DEM  extending  beyond  the  area  of 
interest.  This  may  be  accomplished  in  ArcGIS  as  well.  Note  that  DEM  extents  that  are  too  large 
and/or  resolutions  that  are  too  fine  may  lead  to  computer  memory  issues. 

1.  Open  the  ArcMap  file  thatwasstarted  in  Section  Enori Reference  source  notfound.,  if 
not  already  open. 

2.  Right-click  the  data  frame  ("Downscaled  Soil  Moisture  Estimation")  in  the  table  of 
contents-Mdd  Data...  Browse  for  the  downloaded  DEM  data  and  click  Add.  Do  this 
for  all  the  DEM  grids  that  you  downloaded.  On  the  Tools  toolbar  (shown  below  -  if  it 
not  displayed,  right-click  on  the  grey  space  where  other  toolbars  are  located  and 

select  "Tools"  on  the  drop  down  list),  click  the  world  w  symbol  to  zoom  to  the  full 
extent  of  a  II  the  data. 
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3.  For  convenience,  define  a  group  layer  by  right-click  on  the  data  frame->New  Group 
Layer.  Double-click  the  New  Group  Layer  and  change  its  name  to  Small  DEMs. 
Highlight  all  the  DEMs  again  and  drag  them  into  Small  DEMs  Now,  drag  Small  DEMs 
below  the  coarse  soil  moisture  grid  so  that  the  coarse  soil  moisture  grid  displays 
above  the  DEMs 


E]  &  Downscaled  Soil  Moisture  Estimation 


0  0  CoarseSM_Wet.tif 
□  0  | 


0  0  ASTGTM_N33E0S3_dem.tif 
0  0  ASTGTM_N33E062_dem.tif 
0  0  ASTGTM_N33E061_dem.tif 
0  0  ASTGTM_N32E063_dem.tif 
0  0  ASTGTM_N32E062_dem.tif 
0  0  ASTGTM_N32EQSl_dem.tif 
0  0  ASTGTM_N31E063_dem.tif 
0  0  ASTGTM_N31E062_dem.tif 
0  0  ASTGTM_N31E0Sl_dem.tif 

Display  |  Source  |  Selection  | 


4.  A  histogram  of  the  data  in  the  DEMs  needs  to  be  evaluated.  Therefore,  double-click 
on  each  of  the  DEMs  within  the  table  of  contents  (on  the  "Display"  tab)  to  open  their 
associated  Layer  Properties.  Open  the  Symbology  tab,  and  select  Classified  in  the 
"Show:"  field  to  the  left.  It  will  ask  you  if  you  want  to  compute  a  histogram.  ClickYes. 


5.  Now  that  the  histogram  has  been  computed,  go  back  to  Stretched  in  the  Show  field 
by  clicking  on  it. 
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6.  Right-click  on  the  Color  Ramp  list  box  and  uncheck  Graphic  View.  Then,  select 
Elevation  #1  in  the  drop-down  list.  Click  OK. 

7.  Repeat  this  process  forevery  DEM  thatyou  downloaded. 


8.  Once  you  completed  the  process,  you  should  have  a  view  that  looks  like  this: 
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9.  Now,  mosaic  the  adjacent  DEM  pieces  together.  Open  AtcToolbox  (shown  in  the 
lower  left  in  the  figure  above).  Open  Data  Management  Tools->Raster->Raster 
Dataset->Mosaic  To  New  Raster  Fill  the  form  out  similar  to  the  following  figure.  Ensure 
that  the  Raster  data  set  name  hasthe  extension  *.tif,  because  TauDEM  only  processes 
this  file  forma  t  (1 IFF). 


10.  Close  the  output  message  box. 
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11.  In  the  table  of  contents,  right-click  the  new  layer  name  and  click  Properties...  Go  to 
the  Symbology  tab.  Right-click  on  the  Color  Ramp  and  uncheck  Graphic  View.  Then, 
select  the  Elevation  #1  color  ramp  from  the  drop-down  list.  Click  OK. 
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12.  Up  until  now,  the  DEM  data  has  been  in  an  un-pnojected,  Geographical  Coordinate 
System.  Also,  the  coordinate  sydtem  may  very  well  be  undefined.  Therefore,  right-click 
the  data  frame  “ Downscaled  Soil  Moisture  Estimation"  in  the  table  of  contents,  click 
Properties...  On  the  Coordinate  System  tab,  see  if  the  Current  coordinate  system  is 
defined.  If  so,  click  Cancel.  If  not,  go  to  "Select  a  coordinate  system".  Predefmed-> 
Geographic  Coordinate  Systems->Worid  ->WGS1984  ora  different,  more  applicable 
geographical  coordinate  system  to  area  of  interest. 

13.  Foruse  in  DownscaleConcept,  a  projected  coordinate  system  is  required.  Therefore, 
project  the  new  DEM  to  a  coordinate  sydtem  associated  with  the  specific  area  of 
interest.  For  Afghanistan,  the  WGS  1984  UTM  Zbne  41N  is  applicable.  Open  the 
ArcTooibox,  and  click  on  Data  Management  Too ls->Projectio ns  and  Tranformations- 
>Raster->Project  Raster. 

14.  Place  the  DEM  into  the  Input  Raster  field. 
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15.  Select  an  output  file  path  for  the  projected  raster,  and  ensure  that  the  file  has 
extension  *. tif  because  DownscaleConceot  requires  11FF  fijes  as  input  Note  that 

Ta uDEM  requires  WF_  fijes  as  wejj  if  usng  it  for_  topographic  anajysjs  of_  the  DEM 

chajactenstics, 


16.  Click  the 


symbol  next  to  Output  Coordinate  System. 
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17.  Press  Select...  Navigate  to  Projected  Coordinate  System->UWI->WGS  1984->  WGS 
1984  UTM  Zone  41N.prj  (Remember  this  projection  is  applicable  to  Afghanistan,  not 
necessarily  your  a  pplic  a  tion. 
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Add 


Cancel 


18.  Choose  either  the  CUBIC  or  BIUNEAR  option  as  the  Resampling  Technique  because 
these  are  the  most  suitable  forDEM  data. 

19.  Choose  the  desired  cell-size.  It  is  suggested  not  to  make  the  cellsize  larger.  Making  it 
smaller,  though,  may  increase  processing  time  significantly.  The  figure  below  displays 
the  difference  between  an  unprojected  coarse  grid  and  a  projected  finer  grid 
(remember  that  this  finer  grid  is  no  more  accurate  than  the  original  coarse  grid's 
accuracy  and  can  also  increase  both  the  processing  time  and  run  into  memory 
issuesduring  runtime  that  may  result  in  crashing  your  computer)... 


40 


20.  Ensure  thatthe  final  Project  Raster  form  looks  similar  to  the  figure  below  and  pressOK. 


21.  Add  a  new  group  layer,  name  it  Big  DEMs,  and  place  it  under  the  coarse  soil  moisture 
grid.  Drag  the  two  DEMs  (unprojected  and  projected)  into  the  group  layer. 
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22.  Applying  the  same  color  ramp  to  this  raster  image  as  the  un-projected  one,  the 
display  should  look  like  the  below  figure. 
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C.2  Geo-Processing  with  TauDEM 

TauDEM  is  used  to  calculate  slope  and  contributing  area  for  each  cell  within  the  DEM.  To  run 
TauDEM  follow  the  instructions  below... 

1.  Open  ArcToolbox. 

2.  If  you  haven't  added  TauDEM  Tools,  right-click  in  the  white  space  in  ArcToolbox, 
select  Add  Toolbox...  and  browse  to  Taudem  in  the  Program  Files  directory  (if  you 
installed  it  already).  HighlightTauDEM_Toolsand  clickOpen. 

3.  Within  ArcToolbox,  open  TauDEM  Tools->Basic  Grid  Analysis->Pit  Remove.  Place  the 
projected  DEM  asthe  Input  Elevation  Grid. 

a.  Enter  input  number  of  processes  consistent  with  how  many  processors  your 
computer  has  If  you  don’t  know  how  many  processors  your  computer  has, 
open  your  task  manager  and  go  to  the  Processes  tab.  The  number  of 
windows  within  the  CPU  Usage  History  field  is  the  number  of  processors  your 
computer  has.  The  figure  below  shows4  processors 
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b.  Click  OK.  Running  Pit  Remove  should  take  a  couple  minutes.  If  it  is  taking 
longer,  there  might  be  a  memory  issue  and  a  DEM  with  coarser  resolution  may 
need  to  be  processed  instead. 


4.  Then,  in  TauDEM  Tools,  click  on  Basic  Grid  Analysis->D-lnfinity  Flow  Directions,  which 
tends  to  have  better  results  than  D-8  Flow  Directions  and  produces  slope  and  angle 
of  flow  direction  as  an  output.  Place  the  DEM  with  filled  pits  as  the  input  raster. 


5.  Click  OK.  This  process  may  take  an  hour  or  two  depending  on  the  extent  and 
resolution  of  the  DEM.  With  a  DEM  the  size  used  in  this  example  (9  smaller  grids 
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connected),  a  64-bit  quad-core  computer  processing  using  3  processes  at  2.67  GHz 
took  2.5  hours  to  finish  the  task.  The  following  console  should  pop  up  when  running 
TauDEM. 


C:\Wi  nd  ows\system32\cmd  .exe 

- 

This  estimate  is  very  approximate. 

Run  time  is  highly  uncertain  as  it  depends  on  the  complexity  of  the  input  data 
and  speed  and  memory  of  the  computer.  This  estimate  is  based  on  our  testing  on 
a  dual  quad  core  Dell  Xeon  E5405  2.0GHz  PC  with  16GB  RAM. 

All  slopes  evaluated.  19705678  flats  to  resolve. 

Draining  flats  towards  lower  adjacent  terrain 


6.  Now,  contributing  area  can  be  calculated.  Go  to  TauDEM  Tools,  click  on  Basic  Grid 
Analysis->D-lnfinity  Contributing  Area.  Deselect  "Check  for  edge  contamination"  to 
ensure  that  missing  valuesdo  not  appearin  the  area  of  interest.  ClickOK. 
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7.  If  desired,  run  D-Infinity  Contributing  Area  with  "Check  for  edge  contamination" 
selected  in  order  to  see  how  many  missing  values  due  to  edge  effects  run  through 
the  area  of  interest. 


a.  After  running  the  Contributing  Area  calculation,  zoom  into  the  area  of  interest 
with  the  Contributing  Area  displayed  on  top. 

b.  Double-click  on  the  contributing  area  layer  to  open  its  Layer  Properties.  On 
the  Symbology  tab,  select  Classified  or  Stretched  in  the  Show  field  to  the  left. 
In  the  bottom  right  of  the  window  is  "Display  NoData  as".  Click  on  the  drop 
down  box  next  to  this  and  select  a  bright  color  (red  or  pink).  Click  OK. 
Observe  the  contributing  area  grid  on  the  map. 

c.  If  there  are  no  missing  values  through  the  area  of  interest,  every  upstream 
point  has  been  included  in  analysis  If  not,  it's  ok,  DownscaleConcept  uses 
more  of  a  relative  contributing  area  anyways,  but  you  may  consider  retrieving 
more  of  the  basin  extent  in  the  original  DEM. 
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C.3  Produce  Aspect  Grid 

In  order  to  determine  solar  effects  on  soil  moisture  estimation,  the  aspect  (direction  in  which  a 
slope  is  facing)  is  used  as  input  to  DownscaleConcept.  Calculate  aspect  in  ArcMap  by  following 
the  instructions  below: 

1.  In  ArcMap,  click  Tools->Extensions...  Ensure  that  the  Spatial  Analyst  extension  is 
selected. 


2.  In ArcToolbox,  double-clickSpatial AnalystTools->Surface->Aspect. 

3.  Place  the  projected  DEM  (not  the  one  with  filled  pits)  into  the  Input  Raster  fie  Id. 
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4.  Ensure  that  the  suffix  is  ".tif"  because  DownscaleConcept  accepts  only  TIFF  file 
formats  as  inputs  Click  OK. 


5.  Drag  the  aspect  layer  below  the  coarse  soil  moisture  grid.  Output  will  look  something 
like  this 
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C.4  Clipping  Terrain  Characteristics  to  Area  of  Interest 

After  terrain  characteristics  have  been  calculated,  all  associated  rasters  will  be  clipped  to  the  area 
of  interest  for  modeling  purposes.  Follow  these  instructions: 

1.  In  ArcToolbox,  open  Data  ManagementTools->Raster->RasterPnocessing->Clip. 

2.  There  are  three  grids  that  will  be  clipped:  Aspect,  ContribArea,  and  Slope.  Place  one 
of  these  three  into  the  "Input  Raster"  field  on  the  Clip  dialog.  Place  the  Coarse  Soil 
Moisture  grid  into  the  "Output  Extent"  field.  Change  the  extension  on  the  Output 
Raster  Dataset  to  ".tif'to  ensure  the  output  raster  is  in  the  TIFF  format.  Ensure  that  the 
dialog  lookslike  the  one  in  the  figure  below  and  then  clickOK. 


3.  If  the  clip  does  not  reduce  the  size  of  the  Aspect  /  Slope  /  Contributing  Area  grid, 
then  make  sure  that  Coarse  SM  has  a  defined  coordinate  system  (ArcToolbox->Data 
Management  Tools->Projections  and  Transformations->Define  Projection)  and  repeat 
Step  2. 
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4.  Repeat  Steps  1  through  3  foreach  of  the  three  terrain  characteristics  (Aspect,  Slope, 
Contributing  Area)  to  reduce  the  size  of  each  of  the  grids. 
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